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УДК 532.593 
ИССЛЕДОВАНИЕ ДВИЖЕНИЯ СЛОЕВ ЖИДКОСТИ 


В.В. ТРЕПАЧЁВ 
(Донской государственный технический университет) 


Проведено исследование дифракционного движения двух слоев несжимаемой тяжелой жидкости вокруг вер- 
тикальной круговой опорной колонны. Использовано приближение твердой крышки для выделения воздей- 
ствия только внутренней волны. Результаты исследования распространяются на все возможные масштабы 
волнового движения слоев жидкости как натурного наблюдения, так и лабораторного испытания. 
Ключевые слова: слои жидкости, опорная колонна, сила нагрузки, волна внутренняя, диссипация, ди- 
фракция, межфазное натяжение, приближение твердой крышки. 


Введение. Воздействие солнечного излучения на водную поверхность часто приводит к образо- 
ванию летней устойчивой двухслойной стратификации жидкости в океане. Различные виды по- 
верхностного опреснения жидкости также приводят к разделению одной и той же жидкости на 
два слоя разной плотности. На границе вод разной плотности существует межфазное натяжение, 
что подтверждено наблюдениями [1]. В технологиях очистки и переработки ряда жидкостей ис- 
пользуют двухслойные пакеты жидкости, изолированные от внешней среды твердым дном и 
крышкой. Была предложена модель для точного выделения эффектов внутренней волны с помо- 
щью приближения твердой крышки. Эти положения использованы для обоснования применяемой 
теоретической модели. Возникновение неустойчивости Кельвина — Гельмгольца слоев жидкости 
при их течении под углом друг к другу изучено в [2]. Случай дифракции поверхностной волны на 
круговой опорной колонне в однослойной жидкости в [3] рассмотрен без учета влияния диссипа- 
ции энергии и поверхностного натяжения. Таким образом, из полученных результатов теория вы- 
текает в виде частного случая [3]. 

Постановка задачи. Рассмотрим движение двух слоев идеальной несжимаемой тяжелой жидко- 
сти в отсутствии вихрей. Верхний слой жидкости имеет постоянные плотность р, и толщину №, 


он сверху ограничен горизонтальной жесткой пластиной, называемой твердой крышкой. Нижний 
слой жидкости имеет постоянные плотность р, и толщину В, , он снизу ограничен жестким гори- 


зонтальным дном. Стратификация жидких слоев является устойчивой, что соответствует неравен- 
ству р,<р›. При этом более плотный слой жидкости находится внизу. Учитывается наличие дей- 


ствия малых рассевающих энергию массовых сил, которые пропорциональны скорости движения 
жидкости согласно теории, предложенной Рэлеем [4]. 
Считаем, что верхний и нижний слой имеют одинаковый коэффициент диссипации и>0. 


При наличии сил диссипации Рэлея соблюдается теорема Лагранжа о сохранении потенциально- 
сти движения жидкости [4]. Потенциальное движение слоев жидкости описывается с помощью 
двух потенциалов скорости: Ф, для значений 0<2<А,, Ф, при -В, <2 <0, где 2 - вертикаль- 


ная координата. Начало декартовой системы координат 0хуг выбрано на невозмущенной грани- 


це раздела жидкостей. В направлении оси х распространяется внутренняя гармоническая волна, 
имеющая некоторую круговую частоту колебаний ‹ и волновое число К, 


по = ехр[-КАХ - ®ё)], (1) 

где Е -— время, «> 0, / -— мнимая единица, К = К, — №, , К, >0, К, >0. Длина волны равна 2л / К,. 
Плоская волна, имеющая постоянный амплитудный множитель п, , затухает вдоль направле- 

ния своего распространения, К, — коэффициент затухания бегущей волны. На границе разде- 

ла жидкостей учитывается влияние межфазного натяжения с помощью коэффициента меж- 
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фазного натяжения «о. Внутри жидкости расположен жесткий неподвижный вертикальный 
круговой цилиндр, который занимает следующую осесимметричную область: 
Хх? +у?<а?, -В, <2<В. (2) 
Круговой цилиндр представляет собой препятствие на пути внутренней волны (1). Вся бо- 
ковая поверхность цилиндра рассеивает внутреннюю волну. Поле рассеянных волн неизвестно и 
может быть найдено из решения краевой задачи во внешней области по отношению к круговому 
цилиндру (2). Вводим потенциалы скоростей и форму границы раздела слоев С(х,у,Ё): 


Ф, =$,(х,у,2)е^' , Ф, =ф,(х,у,2)е ' б=и(х, уе”. (3) 


Давление в разных слоях жидкости определяем с помощью интегральных соотношений 
линейной теории: 
Р =-р.[(ю+н)фе” + 92], 0<2<А; Р, =-р.[(1ю+ и)ф,е"" + 92], -В <2<0, (4) 
в которых опущены произвольные постоянные и функции времени [3, 4]. 
Постановка краевой задачи для потенциалов скорости принимает вид: 


5 0 0%, | 


ГоЁ 














дх? ду’ 02 ' и. © 
в № -р.<2<0; (6) 
0х? 0у? 02° " ° ' 
д д . 
ее = штоб У), 250; (7) 
| 0 | 0ф 0 
р» | — мо)ф, - 9 т — р: © — Пиф) ф, — 9 о. о 573 0, 2=0; (8) 
м(х, у) = иее””” +1 (Хх, у); (9) 
12 (х, У) > 0, г? = +у? > о; (10) 
06 . 9 . 
о г=а, Рей: а. г=а, -1, <2<0; (11) 
0$, 0ф> 
— =0, 1=И; — =0, 71=-В.. 
92 . 02 ь Иа 


Уравнения Лапласа (5), (6) получены из условий неразрывности для слоев несжимаемой 
жидкости. Кинематическое условие (7) указывает на непрерывность вертикальной скорости на 
границе раздела жидкостей. В (8) представлено динамическое условие на границе раздела слоев 
жидкостей, которое учитывает влияние эффекта межфазного натяжения наличием коэффициента 
о. Межфазное натяжение вызывает скачок нормального давления на границе раздела. Соотно- 
шение (9) отражает структуру рассматриваемого волнового движения слоев жидкости. Первое 
слагаемое имеет известное значение амплитуды п, и представляет собой бегущую волну, кото- 


рая называется падающей внутренней волной. Соответственно второе слагаемое называется рас- 
сеянной волной п›(х,У). Диссипация энергии и кольцевая расходимость волн вызывают согласно 


условию (10) затухание волнового движения границы раздела для слоев жидкости на бесконеч- 
ном расстоянии от начала координат. Граничные условия на боковой поверхности цилиндра в 
(11) указывают на отсутствие нормальной скорости жидких слоев на ней. Постоянная величина а 
равна радиусу вертикального цилиндра, рассеивающего на своей поверхности волну (1). 
Граничные условия на горизонтальных поверхностях в (12) указывают на отсутствие нор- 
мальной скорости жидкости на них. Граничные условия (11) и (12) называются условиями непро- 
ницаемости для поверхностей, ограничивающих область течения идеальной жидкости. 
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Построение решения. Решение для потенциалов скорости краевой задачи (1) - (12) имеет вид 











функций: 
_ ю@К(2-В,)| , (13) 
ф.(х,У,2) К (АА, ) м(х, У); 
_ юбт[К(2 + В,) | (14) 
ф.(х,у,2) К (АВ, ) "(х,У), 


которые удовлетворяют кинематическому условию (7), граничным условиям непроницаемости на 
горизонтальных плоскостях в виде дна и твердой крышки в (12). Функция т(х, у), заданная со- 


отношением (9), содержит неизвестное слагаемое в виде внутренней волны п,(х, у), рассеянной 
на цилиндре. Рассеянная внутренняя волна ",(х, у) удовлетворяет уравнению Гельмгольца 














2 2 (15) 
г т + Ки, = 0, 

Хх у 

а также двум граничным условиям: 
[№ 16 

ОП > не и гла: (16) 
ог ог 

1>(г,0) > 0, г >. (17) 


Декартовы координаты и полярные координаты в формулах (15) - (17) связаны известны- 


ми соотношениями: х =гс0$0, у=г9п0, г=./х*+у?. 

Волновое число, входящее в соотношения (1), (9), (13), (14), (16) и уравнение (15), удов- 
летворяет дисперсионному соотношению внутренней волны на границе раздела двух слоев жид- 
костей 


(©? — Гон[СЕ(КЯ ›) + усЕ(КА,)] = 9К[(1-у) +ВК?]. (18) 
Определяем параметры у, В , которые входят в дисперсионное уравнение (18), 
т в. в-®. (19) 
Р› р29 
Параметр у имеет смысл относительной плотности слоев жидкости. В случае 
0 <у<1 (20) 


движение двухслойной жидкости устойчиво по Кельвину — Гельмгольцу [2]. Действительно, при- 
нимая и=0, находим из (18), что каждому действительному значению волнового числа К соот- 
ветствует два действительных значения частоты колебаний « при выполнении двойного нера- 
венства (20). Полагая и =0, у=0 в уравнении (18), выводим в виде частного случая дисперси- 


онное соотношение для поверхностных волн [3]. 

В рассматриваемой задаче частота колебаний принимает действительное значение, а вол- 
новое число — комплексное число, удовлетворяющее неравенствам, наложенным на его действи- 
тельную и мнимую части в (1). Экспоненту, входящую в условие (16), следует представить в виде 
тригонометрического ряда согласно формуле 8.5114 из [6] 


е-—*“<°° — у "(Г)", (Кг) с0$(176), (21) 


где 7,(7) — функция Бесселя первого рода [6], =, =1, &, =2, в, =2,.... 

Решение уравнения Гельмгольца (15) отыскиваем в виде разложения по функциям Бессе- 
ля третьего рода (функциям Ханкеля) Н“?(2) [6], которые удовлетворяют условию убывания на 
бесконечности (17) при комплексном значении волнового числа 
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п (г, 6) = У В.Н) соз(пб), г>а. (22) 
п=0 


Неизвестные коэффициенты В, определяются из граничного условия на поверхности ци- 
линдра (16) 


=, (-1)" 7, (8) (23) 


В = Е = Ка, 


, Ц 
: а: 
где штрих означает первую производную по безразмерной переменной & = Ка, которая имеет 
смысл безразмерного волнового числа. Построенное решение уравнения Гельмгольца (15) 
удовлетворяет всем поставленным условиям, т. е. (16) и (17). 
Подставляя (22) и (23) в (9), находим формулу для расчета формы поверхности для гра- 
ницы раздела слоев жидкости (г, 0) в виде ряда в полярной системе координат 








(24) 





м(г, 0) = пе, (И - „(Кг) - о на) веним) „раз 


Формула (24) описывает суперпозицию падающей волны и рассеянной волны. Подставляя 
формулу (24) в соотношения (13) и (14), определяем формулы расчета потенциалов скоростей, 
преобразованные формулы (13) и (14) в каждое интегральное соотношение линейной теории (4), 
находим законы распределения давления в каждом из слоев жидкости. 

Полагая, что г =а в (4), выводим законы распределения давления слоев жидкости на по- 
верхности вертикальной круговой колонны 








= 2 ; с [2 и в, | ГоЁ В 
Р(а, 6,2) =-р, С Пл) аи) п(а,6)е”“" + 92|, 0<2<В; (25) 
СН К(2+1,) а 
Веру в ле) | ] (а, 6де ‘ое — В, <2<0. 


Давление состоит из динамической и гидростатической составляющих. Функция 
"(а,0)е”' в (25) определяет форму возвышения границы раздела непосредственно на самой по- 
верхности цилиндрической колонны 

211. <> „\„ С0$(л6) (26) 
"(а, 0) о = Г)? '(2) ы 
п п=0 Н, (=) 


Формула амплитуды колебаний (26) выводится из (24) с помощью соотношения 8.477, из 








[6]. 

Сила Е , действующая на цилиндрическую опору, находится путем вычисления поверхно- 
стного интеграла от законов распределения давления слоев жидкости на поверхности колонны 
(25) и имеет вид 


2 0 в; 
Е= —а | со] | Р,(а, 0,2) + аби |во7 = 
0 


в о 





— 4(р. — р.) 9% а[(1-у)+ВА*] её 7 = = ка (27) 

&Н, (8) — К[сё®(КЬ,) + усЕИ(КВ, )] | | 
где / -— единичный вектор, направленный вдоль оси х. Решение для случая однослойной жидко- 
сти, которое описывает только воздействие поверхностных волн на вертикальную колонну 
[3], вытекает из полученного решения в виде частного случая, соответствующего значениям 
плотности верхнего слоя р, =0, коэффициента диссипации энергии и =0, коэффициента 


межфазного натяжения а =0. 
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Заключение. 1. Давление слоев жидкости на поверхность колонны (25) испытывает скачок пер- 
вого рода на границе раздела. 

2. Увеличение межфазного натяжения при фиксированном значении волнового числа уве- 
личивает величину модуля силовой нагрузки (27). 

3. Уменьшение скачка плотности р, -р, при фиксированном значении волнового числа 


уменьшает величину модуля силовой нагрузки (27). 
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